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1 .  Introduction 

The  problem  of  electromagnetic  coupling  from  one  region  to  another 
region  through  an  aperture  in  a  ground  plane  can  be  formulated  in  terms 
of  two  admittance  operators,  one  for  each  region  [1].  These  admittance 
operators  are  complex  and  symmetric*  Recently,  a  theory  for  the  charac¬ 
teristic  modes  for  apertures  has  been  proposed  [2]*  These  modes  are  de¬ 
fined  as  the  eigenfunctions  of  a  certain  generalized  eigenvalue  equation 
involving  the  admittance  operators  and  their  real  parts*  Because  of  the 
particular  choice  of  the  eigenvalxie  equation,  the  characteristic  modes 
are  real  (or  equi-phasal)  and  orthogonal  with  respect  to  both  the  admit¬ 
tance  operators,  their  real  parts,  and  their  imaginary  parts  over  the 
slot*  Furthermore,  the  characteristic  fields  produced  by  the  character¬ 
istic  modes  are  orthogonal  over  the  radiation  sphere.  For  small  aper¬ 
tures,  the  characteristic  mode  theory  reduces  to  an  augmented  Bethe  hole 
theory,  i*e.,  the  aperture  is  described  by  a  susceptance  term  related  to 
the  polarizability,  plus  a  conductance  terra. 

In  this  paper,  the  general  theory  is  specialized  to  an  infinitely 
long  slot  in  a  ground  plane  in  an  unbounded  medium.  The  purpose  of  the 
paper  is  to  illustrate  the  theory,  then  apply  it  for  a  complete  analysis 
of  the  slot  problem  vdien  the  slot  is  illuminated  by  a  uniform  TE  (to  the 
slot  axis)  plane  wave.  Specifically,  the  characteristic  currents  and 
fields  are  computed  for  different  slot  widths.  These  are  then  used  to 
compute  the  quantities  and  parameters  of  importance  usually  encountered 
in  electromagnetic  field  compatibility  problems,  such  as  the  equivalent 
magnetic  current  of  the  slot,  transmission  coefficient,  and  transmitted 
field  pattern  both  near  to  and  far  ffom  the  slot.  Analytic  expressions 
for  the  special  case  of  the  narrow  slot  are  also  given. 
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The  infinitely  long  slot  problem  has  been  extensively  considered  by 
researchers  over  the  years.  A  long  list  of  previous  work*  can  be  found 
in  [3]  and  [4].  The  theory  of  characteristic  modes  is  seen  in  this 
paper  to  be  both  general  and  computationally  efficient.  The  extensions 
of  the  work  here  to  a  slot  in  a  ground  plane  separating  contrasting  me¬ 
diums,  and  to  a  uniform  IM  plane  wave  excitation  of  the  slot  are  under 
preparation. 

2,  Basic  Formulation 

Let  the  excitation  be  a  plane  wave  of  unit  amplitude  incident  on  the 
slot  at  an  angle  9  from  the  left  of  the  screen  (see  Figure  1).  This 
wave  is  assumed  to  be  uniform  and  TE  to  the  slot  axis,  and  therefore  has 
the  field  distribution 


-jkCx  sinG  +  z  cos6) 


I 

( 


(1) 


-jicCx  sinG  +  z  cosG) 


where 


K  =  CJ  /(iS 


(2) 


K 


are  the  medim  wave  nmber  and  impedance,  respectively. 


Becaxise  of  the  presence  of  the  slot,  part  of  the  incident  field  is 
reflected  back  into  the  z<0  half-space,  while  the  r^st  of  it  is  trans¬ 
mitted  through  into  the  z>0  half-space.  The  total  field,  incident  plus 
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AN  INFINITELY  LONG  SLOT  IN  A  GROUND  PLANE 
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scattered,  (E,  H) ,  must  have  zero  electric  field  component  tangent  to 
the  screen,  and  continuous  tangential  electric  and  magnetic  fields 
across  the  slot.  A  field  equivalence  theorem  is  used  to  divide  the 
problem  into  two  decoupled  parts. 

Let  the  exciting  field  be  incident  while  the  slot  is  covered  by  a 
perfect  conductor.  This  field,  often  referred  to  as  the  short  circuit 
field  (E^^t  is  given  by 


sc  ^ 

E  =  “2c  |jcos0  sin(icz  cos0)  x 


+  sin0  cos(icz  cos6)  z 


-IKX  sinO 
e 


TfSC  _  r.  f  rs\  sinO 

H  =2  cosCkz  COS0)  e  y* 


(3) 


By  the  field  equivalence  theorem  [5,  Section  3-5],  the  field  to  the  leff 

sc  sc 

of  the  screen  is  identical  with  (E  ,  H  )  plus  the  field  (E (M) ,  H(M)) 
produced  by  the  magnetic  current  sheet 


M  =  z  X  E 


z=0 


(4) 


on  the  slot  idiile  it  is  covered  by  a  perfect  conductor.  The  field  to 
the  right  of  the  screen  is  then  identical  with  the  field  (E (-^M) , 
produced  by  the  magnetic  current  sheet  -M  on  the  slot  while  it  is  cover¬ 
ed  by  a  perfect  conductor.  Figure  2  shows  the  equivalent  situations. 

Since  the  slot  is  tiniform  along  the  y-axis,  and  since  the  incident 
field  has  only  an  component  that  does  not  vary  with  y,  so  does  the 
scattered  field.  Thus,  the  total  field  is  uniform  and  TE  to  the  y-axis, 
and,  as  follows  from  (4),  M  has  only  a  y-component  that  does  not  vary 
with  y: 
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(b) 


(c). 


Original  problem 


Equivalent  situation 
for  t<0 


Equivalent  situation 
for  t>0 


FIG.  2 


THE  EQUIVALENT  SITUATIONS 
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M  =  M(x)  (5) 

The  field  produced  by  M  to  the  left  of  the  screen  is  then  given  by  [5, 
Section  3-12] 


E(M)  =  -  V  X  F(M) 

f 

H(M)  =  -  j  -  F(M) 

) 


(6) 


v^ere  jF(M)  is  the  electric  vector  potential  produced  by  M  in  the  pres¬ 
ence  of  the  complete  screen  [5,  Section  S-”?]: 

w 

r  1  r  (n)  _ 

F(M)  =  |~  I  M(x*)  H  (k/z^  +  (x-x*)^)  dx*  I  y*  (7) 

-  j  0  )  ^ 

“W 


Here  Hq  is  the  Hankel  function  of  the  second  kind  and  zero  order.  Be¬ 
cause  of  the  linearity  of  the  electric  vector  potential,  the  field  pro¬ 
duced  by  -M  to  the  right  of  the  screen  is  negative  of  that  in  (6). 

The  total  tangential  electric  field  clearly  vanishes  at  the  conduct¬ 
ing  screen.  Furthermore,  the  tangential  electric  field  is  continuous 
across  the  slot  by  virtue  of  placing  magnetic  current  sheets  of  opposite 
signs  on  the  opposite  sides  of  the  slot.  The  continuity  of  the  tangen¬ 
tial  magnetic  field  across  the  slot,  however,  requires  that 

w 

^  I  M(x’)  h'^\kU-x'|)  dx'  =  2  xeS  (8) 

C  J  0  '  ' 

“W 

as  readily  follows  from  (3),  (6),  and  (7). 

In  the  next  section,  the  characteristic  currents  of  the  slot  are  de¬ 
fined.  These  are  then  utilized  to  solve  (8)  for  M. 
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3.  The  Characteristic  Currents  of  the  Slot 

The  integral  equation  (8)  can  be  put  in  the  operator  form 

Y(M)  =  I 


where 


w 

Y(M)  =  “"  I  M(x*)  H  (kIx-x*!)  dx* 
^  J  0 

“W 


I  =  2  e’j"* 


•> 


xeS 


(9) 


(10) 


Since  [5,  Appendix  D] 

(j\ 

Hq  ^(O  =  JqCO  -  jNo(0  (11) 

where  Jq  and  are  the  Bessel  functions  of  the  first  kind  and  zero  or¬ 
der  and  of  the  second  kind  and  zero  order,  respectively,  Y(M)  qan  be 
written  as 

Y(M)  =  G(M)  +  jB(M)  (12) 


where 


w 

G(M)  =  -  I  M(x')  J-(k|x-x'|)  dx' 

5  j  0 

-w 


B(M)  = 


w 

M(x*)  Nq(k1x-x*|)  dx*  . 
~w 


f  xeS 


(13) 


Following  Harrington  and  Mautz  [2],  the  characteristic  cxirrents  of 
the  slot  are  defined  to  be  the  eigenfunctions  of  the  eigenvalue  oper¬ 
ator  equation 
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Y(M  )  =  y  G(M  )  (14) 

n  n  n 

and  are  so  normalized  that 

<M  ,  G(M  )>  =  1.  (15) 

n  n 

In  (15),  denotes  the  inner  product 

w 

<C,  D>  =  [  C  (x)  D(x)  dx  (16) 

J 

-w 

* 

\diere  C  (x)  is  the  complex  conjugate  of  C(x).  Put 

y  =  1  +  jb  •  (17) 

n  n 

Then,  using  (12)  and  (17),  (14)  becomes 

B(M  )  =  b  G(M  ).  (18) 

n  n  n 

Clearly,  Y  is  symmetric  (with  respect  to  primed  and  unprimed  coordi¬ 
nates),  and  so  are  G  and  B.  G  and  B  are  also  self-adjoint,  since  for 
any  and 

<Mj^,  G(M2)>  =  <M^,  j  (Y  +  Y*)(M2)> 

=  (Y  +  Y*)(M^) , 

=  <G(M^),  M2>  (19) 

and  similarly  for  B,  Furthermore,  G  is  positive  definite,  since  the 
time- average  power  radiated  into  the  z>0  half-space 
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rad 


~  Re 


w 


E  X  H 


rw 


•  z  dx 


z=0 


Re 


<Y(M),  M> 


<G(M),  M> 


(20) 


is  always  positive.  Tlie  second  equality  in  (20)  follows  from  (4),  (5), 

(6),  (7),  and  (10).  It  is  then  a  standard  practice  [6,  Section  1-25]  to 

prove  that  all  b  -  and  hence  M  ,  are  real  ,  and  that  M  can  be  chosen  to 
n  n  n 

satisfy  the  orthogonality  relationships 


1.  <M  ,  G(M  )> 
ra  n 


2.  <M  ,  B(M  )> 
m  n 


3.  <M  ,  Y(M  )> 
in  n 


nra 


b  6 
n  nm 


(21) 


(1  +  jb  )  6 

n  nm 


where  6  is  the  Kronecker  delta  function  (0  if  in^n,  and  1  if  m=n)  . 
nm 

All  the  currents  on  a  slot  in  an  infinite  screen  in  an  xmbounded  me- 
diuin  are  required  to  radiate  some  power  however  small.  As  can  be  seen 
from  (18)  and  (20),  the  characteristic  currents  corresponding  to  very 
large  b*s  are  basically  non~radiating .  In  the  limit  as  b->^,  these  cur¬ 
rents  cannot  therefore  possibly  exist.  Only  when  the  z>0  region  is 
bounded  by  a  perfect  conductor,  so  that  no  rAdiation  occurs,  these  cur¬ 
rents  are  physically  present.  (G  is  then  positive  semi-definite  rather 
than  positive  definite.)  When  the  slot  is  very  narrow,  this  is  almost 
true,  and  such  currents  can  be  evaluated  under  narrow  slot  approximat¬ 
ions,  as  is  shown  in  a  later  section.  In  any  case,  however,  all  the 


currents  have  to  exhibit  the  edge  property 


[7,  Sec¬ 


tion  1-4]. 


Ol  ■  -  as  x-^±w 

.  [/w*-x*j 


4.  Modal  Solution  of  the  Operator  Equation 

A  modal  solution  of  (9)  for  the  magnetic  current  M  over  the  slot  is 
obtained  in  this  section. 

Put 

M  =  Z  V  M  (22) 

n  n 
n 


\(diere  M^  are  the  characteristic  currents  of  the  slot,  and  are  complex 
coefficients  to  be  determined.  Substituting  (22)  into  (9),  there  then 


results 


Z  V  Y(M  )  =  I.  (23) 

n  ra 
n 


Taking  the  inner  product  of  (23)  with  each  M  ,  one  obtains 

m 


z  V  <M  ,  Y(M  )>  =  <M  ,  I>.  (24) 

n  m  n  m 

n 


Because  of  the  third  orthogonality  relationship  of  (21),  only  the  mth 
term  in  the  summation  survives.  Thus 


V 

m 


<M  ,  I> 
m 

1  +  jb 

m 


(25) 


Substituting  (25)  into  (22),  M  then  becomes 
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<M  ,  I> 

M  =  £  M  (26) 

The  magnetic  current  M  given  by  (26)  is  called  the  modal  solution  of 

(9). 


5.  Power  Considerations 

The  total  complex  power  entering  the  slot  is  basically 


w 

f  * 

P.  =  |ExhI  •zdx. 
in  j  ~  “ 

“W  z=0 


(27) 


Using  (3),  (4),  (5),  (6),  (7),  (9)  and  (10),  (27)  becomes 

P.  =  <1,  M>  -  4  <Y(M) ,  M> 
in  z 

=  j  <1,  M> 

=  j  <Y(M) ,  M>. 


Furthermore,  using  (26),  (28)  becomes 


P. 

in 


-  jb 
1  +  b 


|<I,  M  >1^ 
n 

) 


(28) 


(29) 


A  comparison  of  (20)  and  (28)  then  shows  that  the  complex  power  trans¬ 
mitted  through  the  slot  is  also  given  by  (29). 

A  parameter  sometimes  used  to  express  the  transmission  character¬ 
istics  of  the  slot  is  the  transmission  coefficient  T.  By  definition, 
the  transmission  coefficient  of  the  slot  is  the  ratio  of  the  time- 
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average  power  transmitted  through  the  slot  to  that  incident  on  the  slot 
[5,  Section  7-12],  Using  (1)  and  (29),  T  is  readily  found  as 


T 


1 

4wc  COS0 


I 

n 


1 

1  +  b" 
n 


|<I, 


(30) 


6 •  The  Characteristic  Fields  of  the  Slot 

The  fields  (E  ,  H  )  produced  by  the  characteristic  currents  M  are 
“n  “n  n 

called  the  characteristic  fields  of  the  slot.  Orthogonality  relation¬ 
ships  for  the  characteristic  fields  over  the  radiation  cylinder  can  be 
obtained  from  those  for  the  characteristic  currents  by  means  of  the  com^ 
plex  Poynting  theorem  [5,  Section  1-10].  These  relationships  are  dual 
to  those  for  the  characteristic  fields  of  a  conducting  body,  and  can  be 
derived  in  a  similar  manner  [8].  Thus 


?  j 

Cy 

^  n 


* 

E 

-m 


•k 

H 

-m 


E  dr  =  6 
~n  nm 


H  dT  = 
— 


nm 


Cy 

•*  G 


f  * 
3.  2u  (n  H 
J 

V 


H  -  e  E 
— n 


E  )  dv  = 


b  6 
n  nm 


(31) 


in  duality  with  those  for  a  conducting  body.  In  (31),  Cy^  is  the  radi¬ 
ation  cylinder,  and  the  integration  in  (31.3)  is  over  the  idiole  space. 
Figure  3  shows  the  integration  domain  for  (31). 

The  third  orthogonality  relationship  of  (31)  states  that  the  differ¬ 
ence  between  the  magnetic  and  electric  energy  stored  in  any  characteris- 


THE  INTEGRATION  DOMAIN  OF  (31) 
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tic  field  is  -b/2u  Joules  for  every  one  Watt  of  radiated  power.  The 
characteristic  fields  corresponding  to  positive  (negative)  b  have  pre¬ 
dominantly  stored  electric  (magnetic)  energy,  and  are  therefore  referred 
to  as  the  capacitive  (inductive)  modes  of  the  slot.  The  characteristic 
fields  corresponding  to  b=0  are  called  the  resonant  modes  of  the  slot. 

A  modal  expansion  of  the  field  radiated  by  M  in  terms  of  the  charac¬ 
teristic  fields  E  and  H  can  readily  be  obtained  using  the  modal  expan- 

—n 

sion  (22)  of  M: 


E(M)  =  Z 

V 

E(M  ) 

=  Z 

V 

E 

n 

n 

- n 

n 

n 

— n 

H(M)  =  Z 

V 

H(M  ) 

=  I 

V 

H 

n 

n 

- n 

n 

n 

—n 

(32) 


In  the  Appendix,  it  is  shown  that  the  field  (E (M) ,  H(M))  given  by  (32) 
converges  in  a  least  squares  sense  on  the  radiation  cylinder. 


7 .  The  Narrow  Slot 

An  important  case  that  requires  a  special  consideration  is  that  of 
the  narrow  slot  (2wic<<l)« 

Since 

k|x-x*  I  _<  2wk  <<  1  (33) 

the  Bessel  functions  Jq  and  can  be  replaced  by  their  small  argument 
approximations  [5,  Appendix  D] : 


15 


Jq(k|x-x'  I)  ==  1  - 


.|2 


Nq(k|x-x'|)  =  J  log  ^|x-x' 


(34) 


\diere  log  denotes  the  natural  logarithm,  and  Y=1.7810724,  G(M)  and  B(M) 
are  then  given  by 


G(M)  - 


K 


W 

1  M( X* )  1 1 

^  J  I 

-w 


K  2^ 

^ix-x*rj  dx* 


B(M) 


w 

la  I 

J 

“W 


M(x*  )  log  l-^lx-x*  II  dx* 
k  ^  J 


xeS 


(35) 


Furthermore,  I  does  not  vary  appreciably  over  S,  and  cau  therefore  be 
approximated  by  the  first  two  terras  of  its  Taylor  expansion  about  x=0. 
That  is 

1^2-  j2K  X  sin0,  xeS.  (36) 

Using  (3),  I  then  becomes 

I  +  j-  E^'^CO)  X,  xeS  (37) 

y  Z  z 

where  (0)  is  written  for  the  point  (x,z)=(0,0). 

The  characteristic  currents  of  the  narrow  slot  can  be  determined  by 
solving  (18)  with  G(M)  and  B(M)  given  by  (35),  viz,: 
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w 

-  —  I  M  (x*  )  log  l^lx-x*  II  dx’ 
IT  j  n  \,  I  ) 

-w 


=  b  [m(x*)  |1-  jx-x*  1^]  dx’,  xeS  (38) 
n  j  n  I  ^  ; 

-w 

An  exact  solution  of  (38)  is  possible:  Using  the  identities 


(39) 


it  readily  follows  that 


M 


1 


b 


1 


I 

r 


.  4  J 


(40) 


is  a  solution  pair  of  (38).  Similarly,  it  follows  from  the  identities 


w 

. 

-w 


w 


-W 


,.2 


dx'  = 


I  2 


-IT  X 


(41) 


that 
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U  X 


^2  2  2 
TTK  W 


(42) 


is  another  solution  pair  of  (38).  In  fact,  these  solution  pairs  are  the 
only  possible  solutions  to  (38),  as  can  readily  be  verified  from  the 
identity  [9,  Appendix] 


w 


T  (— ) 
n  w 


log^ 


x-x  I  I  dx 


-ir  log(^)  Tq(^)  if  n=0 


(43) 


-w 


-  -  T  (-) 
n  n  w 


if  n>0 


where  is  the  Chebyshev  polynomial  of  the  first  Hind  and  nth  order, 
and  the  two  identities  in  (39)  and  (41)  involving  logarithms  are  the  n=0 
and  n=l  cases  of  (43),  respectively.  In  (40)  and  (42),  and  are 
constants  to  be  determined  according  to  (15).  Thus 


M 


1 


(44) 


=  -L_  _ x  - 

TTKW^  /W^-X^ 


(45) 


Furthermore,  it  readily  follows  from  (40)  and  (42)  that  the  chare  ter is- 
tic  values  of  the  narrow  slot  are  very  large  positive  numbers  vdiose  ra¬ 
tio  satisfies 


log 


-4 

YKw' 

4 


2  2 
K  W 


4 

YKWv 

4  ^ 


2  2 
K  W 


r2  .2 


(46) 


The  equivalent  magnetic  current  M  of  the  narrow  slot  is  now  given  by 
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M 


<M  ,  I>  <M  ,  I> 

- ± - ; -  M  +  - = -  M 

1  +  jb^  1  1+1^2  ^ 


<M  ,  I> 

— = -  M  -  i 

1  +  jb^  1  J 


<M2.  i> 


y 


V  +  log 


KW  ,  .IT 

rj  *  J  2 


E^^'CO)  X 
z 


(47) 


as  readily  follows  by  substituting  (37),  (40),  (42),  (44)  and  (45)  into 
(26).  In  (47),  v=log(Y)  is  the  Euler’s  constant.  Higher  order  solu- 
tions  can  be  obtained  by  retaining  more  terms  in  the  small  argument  ap¬ 
proximations  of  Bessel  functions  in  (34),  but,  in  view  of  (46),  the  con¬ 
tribution  of  the  higher  order  characteristic  currents  is  negligible. 
Incidentally,  the  magnetic  current  in  (47)  is  identical  with  the  solu¬ 
tion  given  in  [9]. 


® •  Solution  of  the  Eigenvalue  Equation 

An  exact  solution  of  the  eigenvalue  equation  (18)  for  the  char ac ter’* 
istic  currents  is  rather  difficult,  if  at  all  possible.  An  approximate 
solution  has  then  to  be  sought. 

Put 


N 

M  Z  U  P,  .  (49) 

n  -  ,  nk  k 

k=l 


In  (48),  S  are  non- over  lapping  intervals  such  that 
tc 
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N 

2w  =  I  |S  I  (50) 

k=l  ^ 


where  |S  |  is  the  length  of  the  kth  interval  (see  Figure  4).  In  (49), 

each  is  a  real  fixnction,  yet  to  be  specified,  that  vanishes  on 

and  U  ,  are  real  coefficients  to  be  determined.  Substituting  (49)  into 
nk 

(18) ,  it  becomes 


N 

Z 

k=l 


nk 


B(P,  )  =  b 
k  n 


N 

I 

k=l 


U  ,  G(P,  )  +  R 
nk  k 


(51) 


where  R  is  a  residual  term. 

A  Galerkin  solution  [10,  Section  1-3]  of  (13)  can  be  obtained  by 
requiring  that  R  be  orthogonal  to  all  P^,  viz.: 

<P^,  R>  =  0,  =  li2,  ..,N.  (52) 

Thus,  taking  the  inner  product  of  (51)  with  each  P^,  and  enforcing  the 
Galerkin  condition  (52),  there  then  results 
N  N 

Z  U,  <P  .  B(P,  )>  =  b  Z  U,  <P  .  G(P,)>,  fi,=l,2,..,N.  (53) 

,  ,  nk  k  n  .  T  nk  £  k 

k=l  k-1 

In  matrix  form,  (53)  becomes 


B  U  =  b  G  U 
n  n  n 


(54) 


where  G  and  B  are  the  N  by  N  matrices 


G  =  =  [<?^,  G(P^)>] 


(55) 


and  U  is 
n 


the  N  by  1  vector 
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X 


FIG.  4 

THE  PARTITION  OF  THE  SLOT 
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-► 

u 

n 


[U 


nk 


The  constraint  equation  (15)  now  becomes 

-►T  =  -► 

U  G  U  =1 
n  n 


(56) 


(57) 


where  the  superscript  T  denotes  vector  transpose. 

The  solution  of  (54)  subject  to  (57)  determines  in  a  Galerkin  sense 
the  first  N  characteristic  currents  of  the  slot. 


9 •  Evaluation  of  the  Matrices  G  and  B 

The  evaluation  of  the  matrices  G  and  B  constitutes  a  large  portion 
of  the  work  involved  in  the  solution.  An  efficient  evaluation  of  these 
matrices  is  therefore  necessary  for  the  success  of  the  solution. 


The  ilkth  elements  of  G  and  B  are  given  by 


X  X. 

a  k 

f  f 

G..  =  I  Po(x)  dx  I  P  (x*)  J^(k|x-x*|)  dx’ 

j  it  j  R  U 


itk 


‘it-1 


^k-1 


(58) 


ilk 


i  k 

f  f 

J  P^(x)  dx  j  P^(x*)  Nq(k|x-'X*|)  dx* 


where  P  are  so  far  unspecified.  A  particularly  simple  choice  for  P.  is 


f  1  on  S, 


k 


(59) 


0  on  S 


22 


which  corresponds  to  a  pulse  expansion  of  the  characteristic  currents. 
Using  (59),  the  ilkth  elements  of  G  and  B  become 


'ilk 


il  k 

f  dx 
J 

X, 


Jq(k|x*-x*  1  )  dx* 


il-1  k-1 


ilk 


I  dx  I  Nq(k1x-x*1)  dx*  . 

Vi  \-i 


Put 


I  x-x'  I  )  dx' 

J  u 


I  N  (k  I  x-x'  I )  dx' 
J  u 


xeS  . 
I 


(60) 


(61) 


Then,  by  the  first  mean  value  theorem  for  integration  [11,  Section  7-18] 

there  exists  points  x^  ,  x„  eS  such  that 

■^0  ^0 


G  ,  = 

Is  1 

G  , (x_  ) 

ilk 

*  il ' 

ilk 

1 

r 

(62) 

The  evaluation 

of  and  is  completed  by  integrating  the 

Bessel 

functions  J  and  N  over  S  ,  and,  for  that  purpose,  any  quadrature  rule 
U  U  K 


can  be  used.  Thus 
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“  I  ■3i  ■'O 

1=1 


'ft  •  -  i  “o 

1=1 


r 

(  ^ 

K 

r 

1 

O 

X 

^  Pi 

2 

1 

|s,P 

K 

[ 

0 

^  Pi 

2 

J 

J 

1 

(63) 


In  (63),  Q  is  the  order  of  the  rule,  are  its  coefficients,  deter¬ 

mine  the  location  of  its  abscissas,  and 


'k+i 


+  X,  )  T  2 


(64) 


is  the  midpoint  of  S^. 

When  evalxiating  the  diagonal  elements  of  B  ( il=k)  ,  Nq  offers  a  loga¬ 
rithmic  singular ity  at  x=x*  that  requires  par ticular  attention •  Put 

NQp(K|x-x'  I)  =  (Nq  -  Nq^)(k|x-x'  |)  (65) 

vdiere  is  the  singular  part  of  Nq  given  by  its  small  argument  approx¬ 
imation  (34  )•  Then 


I  i 

r  f 

I  dx  I  (N*  +N  )(k|x-x*|)  dx* 

j  j  Os  Op 

X  ,  X  , 

«,-!  l-l 


-  -  |S„M 

TT  ' 


1 

|SJ^ 

1 

2  log 

YK  2 

-  3 

J 

X  X 

1  1 

1  dx  I  N_  (k  I  x-x'  I  )  dx'  . 

J  J  Op 

^-1  ^-1 


(66) 


Nq^  has  no  singularity  at  x=x' ,  and  can  therefore  be  integrated  by  the 
first  mean  value  theorem  and  quadratures.  Ihus 
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|S„  P 

IT  '  ' 


2  log 

IS  n 

'  i ' 

-  3 

CN 

1 

1  J 

1 

J 

- 1 1^1’  -ii  V 
1=1  ^ 


r 

K 

"Nn  ■ 

Op 

1 

i  2 


(67) 


for  some  x.,  eS  . 

N  a 

Op 

Actiially,  finding  such  points  x  ,  ,  and  x^  is  at  least  as  dif- 

^0  ^0  ^Op 

ficult  as  computing  the  integrals  themselves.  For  sufficiently  small 
|S^|,  however,  the  midpoint  of  can  replace  these  points  vdiile  intro¬ 
ducing  nearly  no  error.  Ihtas  ,  eliminating  the  common  factor 


Sk  ■>! 

1=1 


for  all  ,k 


1=1 


(68) 


®Jlk  ~  ' 


-  S, 


Ji 

1 

Is  n 

2  log 

'  i  ' 

-  3 

]Tf 

YK  2 

1 

L 

+  :r  I  q.  N-  (u.  )  / 
2  Op  1  I 


if  «,=k 


where 


f  Is  M 

k 

U.  =  K 

1 

X  ,  - 

(69) 


The  matrices  G  and  B  so  obtained  are  the  same  as  those  would  result  from 
enforcing  the  point  matching  condition 


R(x^^,)  =  0. 


il=l,2,  ..,N 


(70) 
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in  (51)  rather  than  the  Gal er kin  condition  (52),  except  for  a  slight  al¬ 
teration  in  B  ,  as  can  easily  be  established. 


10.  Numerical  Results 

The  characteristic  currents  and  fields,  as  well  as  the  equivalent 
magnetic  current,  radiation  pattern,  and  transmission  coefficient,  have 
been  computed  for  different  slot  widths.  Some  of  the  results  obtained 
for  the  0.4X,  0.5X,  and  l.OX  slots  are  given  in  this  section. 

In  the  actual  computation,  polynomial  expansions  of  the  Bessel  func¬ 
tions  Jq  and  Nq  [12,  Articles  9-4. 1,2,3]  are  utilized,  while  all  the  in¬ 
tegrals  are  computed  using  an  eight-point  Gaussian  quadrature  [12,  Table 
25.4].  In  evaluating  the  far  fields,  the  Hankel  function  is  first  re¬ 
placed  by  its  large  argument  approximation  [5,  Appendix  D]: 


•J? 


(71) 


Then,  using  the  typical  radiation  zone  approximations  [5,  Section  2-10], 


^  =  k/  +  (x-x')^  =  k(p  -  x'cosi|j) 


i  =.  -1 

K  KP 


(72) 


for  all  points  (x,z)eCy  ,  p»2w,  where  is  the  angle  p=xx+zz  makes  with 

■'  p  _  _  _ 

the  x-axis,  the  electric  vector  potential  produced  by  the  characteris¬ 
tic  current  at  any  point  (p,ilj)eCyp  is  readily  found  as 
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F(M  ) 
- n 


k+1 


/2 


JTTKp 


“  j  Kp  N 


k=l 


nk 


j  Kx'  cosilJ 
e  dx* 


1 


N  iCosi> 

J,  "„k  i\i ' 

t  k=l 


Z- 


(73) 


The  far  characteristic  fields  then  follow  by  substituting  (73)  into  (6). 

An  IMSL  Library  2  subroutine  **EIGZF**  [13]  is  used  to  solve  the  ma¬ 
trix  eigenvalue  equation  (54)  for  the  characteristic  values  and  currents 
In  all  the  computer  runs,  a  "performance  index"  has  consistently  come 
less  than  one,  indicating  that  the  subroutine  has  performed  well.  The 
convergence  patterns  for  the  characteristic  values  for  the  0.4x,  0.5X, 
and  l.OX  slots  are  shown  in  Tables  1,  2,  and  3,  respectively.  As  can  be 
seen,  the  convergence  is  always  montone,  either  tipwardly,  or  downwardly. 
Only  for  the  l.OX  slot,  b^  first  decreases  mono  tonic ally ,  but  later  in¬ 
creases,  although  this  can  be  attributed  to  rounding  errors.  Further¬ 
more,  the  convergence  of  the  lower  order  characteristic  values  is  gener¬ 
ally  faster  than  that  of  the  higher  order  ones.  The  convergence  of  the 
dominant  characteristic  current  for  these  slots  is  shown  in  Figures  5, 
11,  and  17. 

The  last  columns  in  Tables  1,  2,  and  3  warrant  an  explanation. 
"EIGZF"  returns  pairs  of  numbers  \diose  quotients  determine  the  charac¬ 
teristic  values.  Both  numbers  are  checked,  and  v^en  the  divisor  of  any 
such  pair  is  found  to  be  an  exact  zero,  an  infinte  value  is  assigned  to 
the  characteristic  value.  Furthermore,  the  power  radiated  by  the  corre¬ 
sponding  characteristic  current  has  always  been  found  to  be  exactly 
zero.  Unfortunately,  it  is  not  true  that  only  a  finite  number  of  char- 
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Table  1.  The  convergence  of  the  characteristic  values  for  a  0.4X  slot. 


N 

"l 

"2 

^5 

•  •  i 

4 

0.3864999 

2.5419659 

52.9851044 

3 

7.2047462x10 

* 

* 

8 

0.3829608 

2.3814417 

39.5235940 

2.7784331x10^ 

4.6240706x10^ 

00 

12 

0.3819787 

2.3403592 

37.0742823 

2.3726822x10^ 

3.3298869x10^ 

00 

16 

0.3815764 

2.3221121 

36.0919691 

2.2324162x10^ 

2.9452827x10^ 

00 

20 

0.3813810 

2.3119735 

35.5731692 

2. 1634724x10^ 

2.7752366x10^ 

00 

->24 

0.3812775 

2.3055932 

35.2564286 

2.1232330x10^ 

2.7000217x10^ 

00 
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Table  2.  The  convergence  of  the  characteristic  values  for  a  0.5X  slot. 


N 

^2 

•  •  i 

4 

0.2486212 

1. 7228569 

22. 1775861 

3 

1.8489054x10 

* 

'k 

8 

0.2514326 

1.6158270 

16.8422341 

0.7267531  xlO^ 

0.7674664x10^ 

00 

12 

0.2524255 

1.5885171 

15.8617172 

3 

0.6228877x10 

0.5547073x10^ 

00 

16 

0.2529714 

1.5763806 

15.4663358 

0.5868004x10^ 

0.4937043x10^ 

00 

20 

0.2533296 

1.5696359 

15.2567972 

3 

0.5690416x10 

0.4661994x10^ 

00 

^24 

0.2535877 

1.5653913 

15.1285194 

0.5586531x10^ 

0.4505765x10^ 

00 

Tab\e  .  The  convergence  of  the  characteristic  values  for  a  1.0  slot. 


29 


• 

• 

• 

■K 

a 

3 

3 

3 

3 

2 

3 

3 

3 

to 

m 

m 

m 

in 

in 

in 

in 

in 

o 

o 

o 

O 

O 

O 

o 

o 

O 

1-H 

r^i 

f~A 

f-H 

f-H 

f-H 

X 

X 

X 

X 

X 

X 

X 

X 

X 

<r 

CM 

m 

00 

<JN 

(3N 

<DN 

o 

CM 

<;*- 

vO 

(50 

LO 

CM 

r— H 

vO 

r— H 

vO 

CM 

r- 

CO 

r-H 

o 

•}c 

VO 

O 

00 

C^ 

CO 

OO 

vO 

CO 

CM 

o 

Mt 

r-* 

CM 

*— f 

in 

VO 

in 

VO 

00 

<r 

CO 

00 

CM 

o 

in 

o 

<r 

o 

(JN 

00 

CO 

• 

• 

. 

. 

. 

• 

. 

• 

. 

in 

■<r 

cn 

CO 

cn 

CM 

CM 

CM 

<r 

<1- 

<r 

<r 

mT 

o 

O 

o 

o 

o 

o 

o 

o 

O 

r-i 

i-H 

•— < 

f-H 

r-H 

f-H 

f-H 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

<!• 

o 

vO 

CO 

o 

VO 

m 

CM 

<r 

m 

VO 

in 

r-H 

•)C 

O 

ON 

CTn 

o 

VO 

in 

JO 

m 

o 

o 

r-. 

o 

CM 

f-H 

vO 

<r 

r-H 

n- 

vO 

m 

00 

ON 

CJN 

m 

O 

00 

VO 

in 

00 

VO 

VO 

VO 

in 

m 

in 

m 

. 

. 

. 

. 

• 

• 

• 

• 

p^H 

o 

o 

o 

o 

o 

o 

o 

o 

cn 

CO 

CO 

m 

CO 

CO 

CO 

CO 

CO 

o 

o 

o 

o 

O 

O 

o 

o 

o 

f— H 

t — 1 

^ — i 

f-H 

r-H 

f-H 

r-H 

X 

X 

X 

X 

X 

X 

X 

X 

X 

<1* 

VO 

■<r 

m 

CO 

<t 

o 

o 

f-H 

o 

Mf 

ON 

r-H 

o 

ON 

VO 

m 

•K 

00 

C7^ 

00 

00 

r-' 

CO 

CO 

vO 

uo 

00 

CO 

CO 

CTN 

VO 

ON 

VO 

VO 

VO 

in 

in 

o 

VO 

CM 

r-H 

00 

t-H 

ON 

00 

00 

r-. 

CM 

CM 

1— > 

1— ^ 

w~~i 

r-H 

r— H 

r-H 

r-H 

• 

• 

■ 

. 

. 

* 

• 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ON 

OV 

in 

CO 

oO 

CO 

in 

VO 

ON 

o 

o 

o 

ON 

CO 

CTN 

f-H 

CM 

CO 

00 

<?■ 

o 

in 

m 

ON 

lOi 

O 

00 

VO 

(ON 

GN 

VO 

r-H 

LO 

LTl 

vO 

(ON 

CM 

m 

1^ 

ON 

CM 

CO 

C7N 

o 

o 

CO 

00 

r-« 

ro 

ON 

<r 

CM 

o 

C3N 

CO 

CO 

• 

• 

• 

• 

• 

. 

• 

• 

• 

• 

m 

CM 

o 

o 

O 

o 

o 

CTN 

CTN 

ON 

CM 

.-H 

F — 1 

r— i 

r-“i 

.—1 

( 

ON 

vO 

00 

in 

f-H 

O 

CO 

cn 

Cv 

r-*-i 

in 

CTN 

<r 

ro 

00 

in 

m 

m 

CM 

VO 

ON 

00 

VO 

o 

VO 

CJN 

00 

00 

<r 

JO 

CM 

ro 

CM 

<t 

m 

(ON 

m 

CM 

o 

CM 

O 

<r 

r-H 

o 

(ON 

00 

CO 

00 

00 

ON 

VO 

m 

in 

in 

<1“ 

<5" 

<T 

Mj’ 

r—f 

i 

1-  H 

*-  i 

r—l 

-H 

f-H 

»— H 

00 

o 

CM 

CO 

<r 

r-H 

ON 

o 

1— ( 

<r 

vO 

vO 

in 

CO 

CTN 

VO 

<r 

.-H 

l-^ 

r- 

VO 

VO 

CM 

CO 

o 

<r 

CM 

O 

o 

o 

o 

o 

m 

CO 

m 

CO 

CO 

CO 

CO 

CO 

CM 

o 

o 

o 

o 

o 

o 

o 

o 

o 

n 

CO 

r^j 

m 

cn 

CO 

CO 

CO 

CO 

o 

o 

O 

o 

o 

C) 

C) 

o 

o 

o 

m 

00 

m 

CM 

r- 

in 

CT> 

<1- 

in 

VO 

o 

c-) 

in 

c> 

O 

CO 

co 

CO 

o 

CM 

VO 

Cs 

CO 

VO 

o 

C^J 

CO 

CM 

<}■ 

00 

ON 

o 

VO 

f-H 

CM 

JO 

ON 

C7N 

m 

o 

in 

CJv 

f-H 

CM 

C^) 

LO 

»— 1 

o 

o 

o 

C) 

O 

r^i 

f-H 

r-H 

o 

o 

o 

o 

C) 

O 

O 

O 

O 

o 

• 

« 

. 

. 

• 

■ 

• 

• 

• 

• 

o 

1 

o 

1 

o 

1 

o 

o 

o 

o 

o 

o 

o 

CM 

VO 

o 

CD 

CM 

VO 

o 

00 

»— 1 

CM 

CM 

CM 

CO 

fO 

30 


acteristic  currents  can  exist  on  the  slot.  Since  single  precision  ar¬ 
ithmetic  is  the  mode  of  operation  used  in  all  computations,  these  char¬ 
acteristic  values  are  actually  0(1/6),  where  6  is  the  smallest  positive 
real  number  different  from  zero  on  the  computer.  Nevertheless,  the  con¬ 
tribution  of  these  currents  is  extremely  small,  and  therefore  haye  not 
been  considered  in  subsequent  computations.  It  appears  that  for  a  slot 
in  a  grotind  plane  in  an  unbounded  medium,  only  a  finite  nmber  of  char¬ 
acteristic  currents  need  to  be  computed.  This  is  also  expected  to  carry 
through  to  slots  in  a  screen  separating  contrasting  raeditims. 

The  computed  characteristic  currents  normalized  to  a  maximum  ampli¬ 
tude  of  unity  and  their  radiation  patterns  for  the  slots  considered  are 
shown  in  Figures  6,  8,  12,  14,  18,  and  20.  This  normalization  is  only 
for  plotting  convenience.  The  equivalent  magnetic  currents  and  radia¬ 
tion  patterns  for  the  slots  are  shown  in  Figures  7,  9,  13,  15,  19,  and 
21,  for  the  case  of  normal  incidence.  It  is  interesting  to  note  that 
the  number  of  lobes  in  each  pattern  is  equal  to  the  order  of  the  charac¬ 
teristic  current  or  field.  Finally,  Table  4  gives  the  ratios  of  the 
powers  radiated  by  the  characteristic  currents  to  that  radiated  by  the 
dominant  characteristic  current  for  each  slot,  whereas  the  transmission 
coefficient  for  the  slots  for  the  case  of  normal  incidence  are  given  in 
Table  5,  all  is  evaluated  with  N  as  specified  in  the  arrow  marked  rows 
in  Tables  1,  2,  and  3.  The  entries  in  Table  4  suggest  that  the  radia¬ 
tion  pattern  for  the  slot  is  bas ically  the  same  as  that  for  its  dominant 
characteristic  current.  This  is  indeed  the  case,  as  Figures  10,  16,  and 


22  readily  show. 
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Table  4.  The  ratio  of  the  power  radiated  by  each  characteristic  current 
to  that  radiated  by  the  dominant  characteristic  current  for  a 
(a)  0.4X  slot,  (b)  0.5X  slot,  and  (c)  l.OX  slot. 

(a) 


n 

P  /  P, 
n  1 

1 

1.0 

2 

0.1595670 

3 

5.1706787x10’^ 

4 

6.8070769x10’^ 

5 

5.3741932x10’^ 

•  •  • 

0.0 

(b) 


n 

P  /  P, 
n  1 

1 

1.0 

2 

0.2678140 

3 

1.1451906x10’^ 

4 

2.3665750x10’^ 

5 

2.9297482x10’^ 

•  •  • 

0.0 
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(c) 


n 


P  /  P 
n 


1  1.0 

2  0.1496020 

3  0.0390829 

4  3.7669604x10' 

5  1.7264361x10' 

6  5.2402006x10 

7  1.1263396x10 

. .  0.0 
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Table  5. 


The  transmission  coeffiecient  for  the  0.4X,  0.5X,  and  l.OA 

slots  for  normal  incidence. 


2w 
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0.4X  0.9757939 

0.5X  0.9777421 


l.OX 


0.9994257 


0.00  0.25  0,.50  0.75  1.00 


Figure  S’  The  convergence  of  the  dominant  chracteristic 

current  for  a  0.4 A  slot  normalized  for  a  maximum 
amplitude  of  unity. 
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Figure  6.  (a)  The  second,...,  (d)  fifth  characteristic 

currents  for  a  0.4  A  slot  normalized  for  a  maximum 
amplitude  of  unity . 
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Figure  7s  The  equivalent  magnetic  current  for  a  0.4X  slot, 
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(a) 

Figure  8:  The  radiation  patterns  of  (a)  the  first, (e) 
fifth  characteristic  fields  for  a  0.4X  slot. 
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Figure  9*.  The  radiation  pattern  for  O.^-X  slot. 
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Figure  10:  Comparison  of  the  radiation  pattern  for  a  0.4 X 
slot  (»  =0.0)  to  that  of  its  dominant  mode. 
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Figure  11 s  The  convergence  of  the  dominant  characteristic 

current  for  a  O.5X  slot  normalized  for  a  maximum 
amplitude  of  unity. 
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Figure  12:  (a)  The  second, (d)  fifth  characteristic  oo 

currents  for  a  0.5^  slot  normalized  for  maximum 
amplitude  of  unity. 
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Figure  13:  The  equivalent  magnetic  current  for  a  0.5X  slot. 

Only  the  (5)  physical  characteristic  currents  are 
used  in  the  summation. 
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Figure  14:  The  radiation  patterns  of  (a)  the  first, (e) 
fifth  characteristic  fields  for  a  0.5^  slot. 
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Figure  15:  The  radiation  pattern  for  a  0.5X  slot. 
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Figure  l6;  Comparison  of  the  radiation  pattern  for  a  0.5X 

slot  (6 =0.0) to  that  of  its  dominant  characteristi- 
c  field. 
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Figure  17 ;  The  convergence  of  the  dominant  chracteristic 

current  for  a  1.0  A  slot  normalized  for  a  maximum 
amplitude  of  unity. 
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Figure  18:  (a)  The  second (f)  the  seventh  characteristic 
currents  for  a  1.0 X  slot  normalized  for  a  maximum 
amplitude  of  unity. 
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Figure  19'  The  equivalent  magnetic  current  for  a  1. slot. 
Only  the  (?)  physical  characteristic  currents  are 
used  in  the  summation. 
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Figure  20:  The  radiation  pattern  of  (a)  the  f irst , . . . , (g) 
seventh  chracteristic  fields  for  1.0 A  slot. 
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Figure  21:  The  radiation  pattern  to  a  l.OX  slot. 
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Figure  22 t  Comparison  of  the  radiation  pattern  for  a  l.OX 

slot  (9=0.0)  to  that  of  its  dominant  chracteristi- 


c  field. 
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11.  Discussion 

The  theory  of  char  ac  ter  is  tic  modes  for  slots  has  been  applied  for 
the  solution  of  the  problem  of  an  infintely  long  slot  in  a  ground  plane 
in  an  unbounded  medium  il linn ina ted  by  a  uniform  TE  (to  the  slot  axis) 
plane  wave.  The  theory  is  a  specialization  of  the  general  theory  of 
characteristic  modes  for  apertures,  and  can  be  applied  to  slots  in  a 
ground  plane  separating  contrasting  mediums  and  illuminated  by  general 
plane  waves • 

The  characteristic  currents  and  fields,  as  well  as  the  equivalent 
magnetic  current,  transmission  coefficient,  and  radiation  pattern,  of 
the  slot  have  been  computed  for  different  slot  widths.  The  theory  is 
rather  general  and  computationally  efficient,  applicable  to  both  narrow 
and  wide  slots.  For  narrow  slots,  the  modal  representation  of  the 
equivalent  magnetic  current  is  readily  recognized  as  an  augmented  multi¬ 
pole  expansion  of  the  cuurent.  Furthermore,  it  appears  only  a  finite 
number  of  characteristic  currents  need  to  be  computed.  The  applicabili¬ 
ty  of  the  theory  is  not  necessarily  confined  to  slots  of  small  or  inter¬ 
mediate  width  compared  to  the  wavelength.  In  fact,  it  is  computationally 
more  attractive  than  a  great  many  other  nimerical  solutions.  This  is 
evident  from  its  performance  in  the  present  situation.  All  the  inter¬ 
esting  features  displayed  here  are  expected  to  carry  through  for  general 
apertures  and  medims. 
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Appendix 


In  this  Appendix,  it  is  shown  that  the  convergence  of  the  field 


E(M)  =  E  V  E  (1) 

-  n 

n 


over  the  radiation  cylinder  is  in  a  least  squares  sense.  The  proof  for 


H(M)  is  similar. 


Put 


N 

E(M)  =  Z  U  E  +  R„(E)  (2) 

-  ,  n  -n  -N  — 

n=l 


where 


N 

IL,(E)  =Z(V-U)E  +  Z  V^E 
n=l  n=N+l 


(3) 


is  a  residxial  term.  Then,  the  convergence  of  E(M)  in  a  least  squares 
sense  over  the  radiation  cylinder  is  equivalent  to  the  requirement  that 
the  norm  of  the  residual 


I  1^  =  I  R^(E)  .  ^(E)  dT 

Cy 


(4) 


be  minimum  if  and  only  if  U^=  n=l,2,,.,N.  Substituting  (3)  into  (4) 

and  using  the  orthogonality  relationship  (31,1),  it  becomes 


|Rj^(E)||  = 


N 

Z 

n=l 


|V  -  U  i 
’  n  n ' 


n=N+l 


(5) 


The  result  then  follows. 
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